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ABSTRACT 


Analytical  and  nximerical  models  are  developed  to  predict 

the  longshore  sand  transport  distributions  across  the  surf 

zone.  The  models,  which  also  predict  the  root  mean  square 

wave  height,  and  the  longshore  current,  V,  are  compared 

with  field  data  acquired  from  Leadbetter  Beach,  Santa  Barbara, 

California,  during  the  intensive  storm  period  in  February,  1980. 

The  breaker  coefficient,  B,and  the  bed  shear  stress  coefficient, 

Cj,  when  equal  to  1.1  cind  0.005,  were  found  to  give  the  best 

agreement  between  the  predicted  H  and  V  with  the  field 

rms 

measurements.  The  model  for  a  plane  sloping  beach  predicts 
the  maximum  sand  transport  at  offshore  distance  X  equal  to 
Q.8X^,  where  Xj^  is  the  mean  breakerline  location,  which  agrees 
with  Komar  (1977a) .  The  longshore  sand  transport  formula, 
suggested  by  CERC  (1977) ,  is  used  to  calibrate  the  model  which 
requires  the  empirical  transport  coefficient  to  be  0.18. 

The  sand  transport  cross-shore  distribution  predicted  by 
the  model  does  not  agree  well  with  the  corresponding  trans¬ 
port  inferred  from  the  field  measurements.  Reasons  offered 
to  explain  the  differences  are  (1)  that  the  model  does  not 
include  the  swash  zone  where  the  maximum  transport  is  found 
and  (2)  that  beach  profile  information  beyond  a  distance  of 
100  m  offshore  was  incomplete.  The  results  indicate  the 
importance  of  including  the  swash  zone  to  describe  the  effect 
of  the  long  waves  in  predicting  the  sediment  transport. 
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I .  INTRODUCTION 


Water  waves  are  one  of  the  principal  causes  of  shoreline 
changes.  When  waves  break  along  the  shore,  they  release  their 
energy  and  momentum  and  give  rise  to  a  longshore  current. 

The  longshore  current,  along  with  the  stirring  action  of  the 
waves,  is  the  primary  mechanism  for  longshore  sediment  trans¬ 
port.  The  longshore  sand  transport  rate  is  an  essential  fac¬ 
tor  determining  erosion  or  accretion  along  a  coast.  If  the 
longshore  sand  transport  rate  can  be  accurately  estimated,  a 
quantitative  picture  of  shoreline  evolution  can  be  evaluated, 
including  changes  of  the  shoreline  due  to  marine  structures. 

In  the  present  study,  analytical  and  numerical  models  are 
developed  based  on  a  longshore  current  model  for  random  waves 
by  Thornton  and  Guza  (1983)  and  a  sediment  transport  formula¬ 
tion  by  Thornton  (1973)  to  predict  the  cross-shore  sediment 
transport  distribution  and  to  compute  the  total  volume  of  sand 
transport  rate.  The  model  is  compared  with  the  field  data 
acquired  from  Leadbetter  Beach,  Santa  Barbara,  California 
( January-February ,  1980) . 

The  experiment  at  Santa  Barbara  was  conducted  as  part  of 
the  Nearshore  Sediment  Transport  Study  (NSTS) .  The  objectives 
of  NSTS  are  to  develop  improved  techniques  and  engineering 
formulae  to  predict  sediment  transport  on  beaches  having 
straight  and  parallel  contours  by  utilizing  field  measurements 
of  waves  and  currents  at  several  nearshore  sites. 


4 


During  the  Santa  Barbara  experiment,  beach  profiles  were 
measured  daily  along  with  extensive  measurements  of  waves  and 
currents.  A  severe  storm  attacked  the  site  of  the  experiment 
during  16-20  Feb.,  1980  causing  significant  changes  in  the 
shoreline  configuration  and  beach  profiles.  Despite  the 
destruction  of  most  of  the  instruments  during  the  storm, 
a  very  substantial  and  valuable  record  of  waves  and  currents 
was  acquired  which  makes  the  Santa  Barbara  data  very  unique  and 
significant  for  nearshore  sand  transport  studies. 

The  present  work  emphasizes  the  data  acquired  on  17  and 
18  February,  1980.  A  strong  longshore  current  was  observed 
due  to  storm  waves  approaching  the  beach  at  moderately  large 
angles  (6®  at  the  breaker  line)  causing  significant  erosion 
in  the  beach  profile.  The  daily  beach  profiles  during  the 
experiment  (January-February)  show  no  well  defined  bars  or 
troughs  and  essentially  have  a  stable  point  95  m  offshore, 
which  did  not  appear  affected  by  erosion  or  accretion.  Thus, 
it  is  hypothesized  that  the  changes  in  the  beach  profile  are 
due  to  sand  transport  in  the  longshore  direction  alone,  and 
the  on-offshore  transport  is  negligible.  Furthermore,  the 
longshore  sand  transport  distribution  across  the  surf  zone, 
predicted  by  the  model,  is  assumed  to  have  a  distribution 
similar  to  the  net  longshore  sand  transport  calculated  from 
the  profile  changes. 
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II.  LONGSHORE  SEDIMENT  TRANSPORT 


A.  GENERAL 

When  waves  approach  the  shoreline  at  an  oblique  angle  and 
break,  they  generate  a  longshore  current  parallel  to  the 
shoreline.  This  current  in  turn  interacts  with  the  waves  to 
produce  a  longshore  sand  transport.  The  movement  of  sand 
particles  is  usually  divided  into  two,  not  always  distinct, 
modes  of  transport:  bed  load  and  suspended  load  transport. 

The  bed  load  is  the  result  of  sand  grains  being  rolled  along 
the  bottom  by  the  shear  of  water  moving  above  the  sediment  bed 
and  is  maintained  by  grain- to- grain  contact.  Bed  load  is 
usually  considered  to  be  confined  to  within  10  cm  of  the  bottom. 
In  the  suspended  load  mode,  the  generally  finer  grains  are 
transported  by  currents  after  they  have  been  lifted  from  the 
bed  by  turbulent  action  of  the  wave-induced  vertical  motion. 

The  total  sand  transport  is  the  sum  of  the  contribution  by 
bed  and  suspended  load  transport.  The  rate  at  which  the  total 
volume  of  sand  is  transported  parallel  to  the  shoreline  is 
termed  "longshore  sand  transport;  Q £" •  The  two  most  success¬ 
ful  and  widely  used  approaches  for  estimating  the  longshore 
sand  transport  are  an  empirical  wave  power  model  and 
Bagnold’s  energetics  model. 

B.  WAVE  POWER  MODEL 

In  the  wave  power  model,  the  rate  of  total  longshore 
transport  is  simply  assumed  proportional  to  the  longshore 


component  of  wave  power  (P^^)  at  the  breaker  line  to  some  power 
n  (Watts,  1953a) . 


where  and  n  are  coefficients  to  be  determined  empirically. 
Assviming  straight  and  parallel  depth  contours,  the  onshore  energy 
flux  in  the  direction  of  wave  propagation  per  unit  length  of 
the  beach  is  EC^  cos  a ,  where  E  is  the  energy  density  of  the 
waves,  Cg  is  the  group  velocity  of  the  incoming  waves,  and 
a  is  the  angle  of  the  wave  crests  relative  to  the  shoreline. 
Then,  the  longshore  component  of  wave  energy  flux  per  unit 
length  of  the  beach  is  given  by: 

=  (ECg) cos  sin  aj^  (2.2) 

where  the  subscript  b  refers  to  the  breaker  line.  Applying 
linear  theory. 


"  F  P  Cg  cos  sin  Ujj  (2.3) 

where  p  is  the  water  density,  g  is  the  gravitational  con¬ 
stant  and  H  is  the  wave  height. 

Watts  (1953a)  obtained  the  first  field  measurements  by 
which  the  sand  transport  rate  could  be  related  to  the  local 
wave  characteristics.  Caldwell  (1956)  combined  additional 
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data  with  those  of  Watts  and  found  the  coefficient  (n)  in 
(2.1)  to  be  less  than  unity,  and  with  a  different  proportion¬ 
ality  coefficient  (K^)  than  Watts. 

Inman  and  Bagnold  (1963)  interpreted  field  and  laboratory 
data  and  obtained  a  linear  relationship  between  and  P^, 
i.e.,  n  equal  to  unity,  which  takes  into  consideration  the 
immersed  weight  of  beach  material. 


_ K _ 

(l-a) (Pg-P) 


(2.4) 


where  is  the  density  of  sand  transport,  a  is  the  porosity 
of  beach  sediment  (taken  as  0.40)  and  K  is  a  dimensionless 
proportionality  coefficient. 

Komar  and  Inman  (1970)  utilized  sand  tracers  to  measure 
longshore  sand  transport  rates.  Using  this  new  data  combined 
with  data  from  Watts  (1953a),  Caldwell  (1956)  and  others,  they 
obtained  0.77  for  the  proportionality  coefficient  K  by  relating 
the  immersed  weight  transport  rate  1^^  to  Pj^ 


I 


I 


0.77  P^ 


(2.5) 


3  3 

For  quartz-density  sand,  =  2.65  x  lo  Kg/m  ,  K  =  0.77, 
equation  (2.4)  becomes, 


=  6.8  P^ 


(2.6) 
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where  in  (2.6)  is  measured  in  units  of  m^/day  and  in 
watts/meter.  Komar  and  Inman  use  the  root  mean  square  wave 
height  in  (2.5)  and  (2.6)  instead  of  the  significant  wave 
height,  suggested  by  CERC  (1977) ,  to  evaluate  the  term. 

The  relationships  described  in  equations  (2.5)  and  (2.6) 
are  almost  purely  empirical  with  no  real  consideration  of 
the  sand  transport  mechanism.  In  addition,  the  coefficient 
0.77  in  equation  (2.5)  has  a  high  degree  of  uncertainty 

Longuet-Higgins  (1972),  in  attempting  to  rationalize  (2.2), 
points  out  that  since  the  energy  flux  (EC^)  is  a  vector  rather 
than  a  second  order  tensor,  the  longshore  component  of  wave 
power  (P^)  can  be  written  equivalently  as  the  product  of  two 
physically  meaningful  quantities; 

C 

Pjj^  =  [E  cos  a  sin  a]C  =  C  (2.7) 


where  is  the  longshore  component  of  the  radiation  stress 

and  C  is  the  phase  speed  of  the  waves.  The  radiation  stress 

(S  )  is  shown  by  Longuet-Higgins  (1972),  among  others,  to  be 
xy 

the  principal  driving  force  for  longshore  currents,  and  thus 

would  be  expected  to  be  important  in  describing  longshore 

sediment  transport.  Vitale  (1981)  compared  the  radiation 

stress  (S^y)  and  the  wave  power  (P^)  as  predictors  of  the 

longshore  sediment  transport  rate  by  conducting  three-dimensional 

movable-bed  laboratory  tests  and  he  found  no  major  difference 

in  choosing  S  over  P.  to  predict  the  longshore  sediment 
xy  jc 

transport . 
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C.  BAGNOLD'S  ENERGETICS  MODEL 


Bagnold  (1963)  relates  the  rate  of  immersed  weight  trans¬ 
port  (Ij^)  to  the  work  done  by  waves  and  current.  The  mechanics 
of  sand  particle  movement  are  described  as  a  back-and-forth 
motion  under  the  action  of  the  wave  orbital  motion  with 
essentially  no  net  transport.  Wave  energy  is  expended  in 
supporting  and  suspending  the  moving  sand  above  the  bed. 

Once  the  sediments  are  suspended,  the  presence  of  a  unidirec¬ 
tional  current  superimposed  on  the  to-and-fro  motion  can  then 
produce  a  net  drift  of  sediment.  Bagnold  (1963)  derived  the 
relationship; 


where  i„  is  sediment  transport  per  unit  width  in  the  direc- 
tion  0  determined  by  the  unidirectional  current  Ug ,  u  is  the 
available  wave  power  supporting  the  sediments  above  the  bot¬ 
tom,  u^  is  the  orbital  velocity  of  wave  motion,  and  K'  is  a 
dimensionless  coefficient. 

Inman  and  Bagnold  (1963)  specified  the  various  parameters 
of  waves  at  the  breaker  line  and  obtained; 


I^  =  K'  (EC  )j^  —  cos 

in 


(2.9) 


where  u  is  the  maximum  horizontal  orbital  velocity  of  the 
m 

waves.  Equation  (2.9)  can  be  viewed  as  a  general  relationship 
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because  it  does  not  specify  the  cause  of  the  longshore  cur¬ 


rent  Vj^;  it  can  be  either  tidal,  currents  of  a  cell  circula¬ 
tion,  wind  generated,  or  due  to  oblique  wave  approach.  Komar 
and  Inman  (1970)  utilized  their  littoral  drift  measurements 
to  test  equation  (2.9)  and  required  K'  to  be  equal  to  0.28 
to  agree  with  their  data.  In  terms  of  and  for  quartz 
density  sand,  equation  (2.9)  is  given  by 

V, 

=  2.5  (EC  )j^  ^  (2.10) 

’  m 

where  is  in  m^/day  and  the  wave  and  current  parameters 
are  in  mKs  units  (watts/meter)  (Komar,  1983) .  Wang  and  Chang 
(1979)  found  good  agreement  between  their  measurements  along 
the  bayshore  of  a  barrier  island  in  the  Gulf  of  Mexico,  and 
equation  (2.9)  for  a  coefficient  K*  =  0.18,  where  the  long¬ 
shore  currents  generally  are  due  to  local  winds  and  tides. 

Bowen  (1981)  applied  Bagnold's  equation  to  the  problem  of 
on-offshore  sediment  transport  on  beaches  and  then  studied 
the  special  case  of  a  normal  incident  wave  and  steady  on- 
offshore  currents  (no  longshore  currents) .  The  waves  are 
described  using  Stokes  second  order  wave  theory  and  the  mean 
on-offshore  cxirrent  using  Longuet-Higgin '  s  (1953)  bottom 
streaming  solutions.  The  model  predicts  an  equilibrium 
beach  profile  as  a  function  of  the  incident  wave  character¬ 
istics.  Bowen's  results  support  observations  that  steep 
beaches  are  generally  coarse  grained,  and  short  period  waves 
cause  a  beach  to  erode. 
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Figxjre  (2.1).  Schematic  of  the  Bagnold  model  for  sand 
transport  wherein  the  orbital  velocity 
uq  due  to  the  waves  places  the  sand  in 
motion,  and  the  current  u^^  provides  a 
net  transport  of  sand,  1^. 
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Ballard  and  Inmem  (1981)  independently  derived  a  very 
similar  model  to  Bowen's  (1981)  to  describe  the  instantaneous 
transport  in  the  nearshore  zone,  including  both  on-offshore 
movements  as  well  as  longshore  movement.  An  important  part 
of  the  study  is  the  effect  of  local  bottom  slope  on  sand 
transport  rate,  as  this  permits  the  analysis  of  the  formation 
of  the  longshore  bars  and  other  inshore  topography.  Bailard 
(1982) ,  extended  the  model  to  predict  the  total  load  trans¬ 
port  for  time-varying  flow  over  an  arbitrarily  sloping  planar 
bed.  The  model  predicts  the  local  equilibrium  beach  slope 
as  well  as  the  local  near  bottom  sediment  transport  rate  as 
a  function  of  the  near  bottom  velocity.  Surf  zone  conditions 
are  described  as  steady  longshore  and  on-offshore  currents 
in  addition  to  a  wave  induced  current  having  a  local  wave 
angle . 

A  number  of  other  models  have  been  proposed.  For  instance, 
Madsen  and  Grant  (1976)  adapt  the  Einstein-Brown  sediment 
transport  equations  to  time-varying  sand  movements  under  com¬ 
bined  waves  and  unidirectional  currents.  The  bottom  stress 
(x)  under  the  combined  oscillatory  motions  of  the  waves  and 
any  superimposed  currents,  and  the  resulting  sand  transport, 
vary  with  time.  The  model  is  relatively  complicated  and  re¬ 
quires  pareuneters  for  which  no  data  are  available.  The 
success  of  the  wave  power  formulation  and  Bagnold's  energetics 
model  is  due  to  their  relative  simplicity  requiring  only  one 
or  two  parameters  specified  from  the  data. 
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In  general,  the  predictions  of  the  littoral  drift  have  a 
high  degree  of  uncertainty.  This  partly  results  from  our 
basic  inability  to  make  accurate  measurements  of  sand  trans¬ 
port  on  beaches.  However,  the  two  approaches  of  equations 
(2.5)  and  (2.9)  can  be  used  to  give  at  least  a  qualitative 
description  of  littoral  drift. 

D.  DISTRIBUTION  OF  LONGSHORE  TRANSPORT  ACROSS  THE  SURF  ZONE 
Bagnold's  (1963)  approach  has  been  also  utilized  to  pre¬ 
dict  the  distribution  of  longshore  sand  transport.  Thornton 
(1973)  analyzed  sand  transport  inside  and  outside  the  surf 
zone  where  the  volume  transport  rate  per  unit  width  (q^)  is 
given  by: 


B, 


.(1  -  ^) 


1/2 


u. 


m 


3  (ECg) 

9x 


(2.11) 


where  B^  is  an  empirical  dimensionless  coefficient.  Thornton 

(1973)  obtained  field  data  at  Fernandina  Beach,  Florida,  to 

evaluate  B  in  equation  (2.11).  The  transport  rate  (q  )  was 
s  s 

measured  using  a  series  of  bed  load  traps  placed  in  a  line 
across  the  nearshore.  The  results  indicate  that  the  maximum 
longshore  transport  is  in  the  breaker  zone  with  B^  =  0.08 
and  stronger  transport  occurs  over  the  bars  rather  than  in 
the  troughs  (Fig.  2.2). 

Komar  (1976b, 1977a)  similarly  utilized  Bagnold’s  (1963) 
model  to  analyze  the  distribution  of  the  littoral  drift  in 
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F 


DISTANCE  FROM  BASELINE  (ft) 


Fi^UT'e  (2.2)  .  The  distribution  of  the  longshore  send 
transport  measured  by  Thornton  (1973) 
at  Fernandina  Beach,  Florida,  by  bed 
load  traps  in  the  positions  as  shown. 
The  dashed  lines  are  based  on  equation 
2.11  with  B  =  0.08.  [from  Thornton 
(1973)]  ® 


23 


terms  of  The  local  immersed-weight  of  sand  transport 

rate  per  unit  of  width  inside  the  surf  zone  (i  )  is  given 

X 

by; 


i(x)  =  (0.5  Cj)  p  h(x)  v(x)  (2.12) 

where  is  the  drag  coefficient  of  oscillatory  wave  motions 
and  h  and  v  are  the  local  water  depth  and  longshore  current,  respec¬ 
tively;  both  are  functions  of  the  distance  x  offshore  from  the  shore¬ 
line.  is  a  dimensionless  proportionality  factor  deter¬ 

mined  by  calibrating  the  model  equation  (2.12)  to  yield 
calculated  from  equation  (2.5)  where  the  total  immersed- 
\  jight  transport  can  be  evaluated  by  integrating  i(x)  across 
the  svirf  zone; 


I„  =  /  i(x)  dx 

^  0 

where  x^  is  the  surf  zone  width  as  shown  in  Fig.  (2.3). 

Komar  (1977a)  calculated  the  distribution  of  both  the  long¬ 
shore  current  and  the  longshore  sand  transport,  Fig.  (2.4), 
based  on  the  longshore  current  analysis  by  Longuet-Higgins 
(1970) .  The  distribution  shows  the  maximum  longshore  current 
at  approximately  (x/Xj^  =  0.65)  and  the  maxim\im  longshore  sand 
transport  at  approximately  (x/Xj^  =  0.80). 

Sawaragi  and  Deguchi  (1979)  utilized  circular  bed  load 
traps  which  enable  them  to  measure  the  on-offshore  sand 
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movements  as  well  as  the  longshore  transport  distribution 
across  the  surf  zone.  The  results  showed  sediment  transport 
to  be  relatively  independent  of  wave  steepness  but 

dependent  on  the  sand  grain  size.  The  maximum  sand  trans~ 
port  is  shown  (Fig.  2.5)  to  occur  at  approximately 
(x/Xj^  =  0.6)  and  maximum  longshore  current  at  (x/Xj^  =  0.4)  . 

Bailard  (1982)  extended  his  earlier  wor)c  in  which  he  used 
second  order  Stokes  theory  to  describe  the  wave  motion.  He 
included  the  Ostenddorf  and  Madsen  (1979)  longshore  current 
model  to  predict  the  distribution  of  sediment  transport  rate 
for  different  ratios  of  the  velocity  of  the  flow  to  the  sedi¬ 
ment  fall  velocity  (Fig.  2.6).  He  finds  the  maximum  long¬ 
shore  current  at  approximately  (x/Xj^  =  0.6)  with  the  maximum 
longshore  sand  transport  at  (x/Xj^  *  0.9). 

Results  from  different  studies  examining  the  longshore 
current  and  longshore  sand  transport  distributions  are  given 
in  Table  I . 
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Figure  (2.5).  Laboratory  measurements  of  the  distribu¬ 
tion  of  the  longshore  current  and  the  sand 
transport  rate  showing  the  positions  of 
the  maximum  sand  transport  and  maximum 
longshore  current  relative  to  the  breaker 
zone.  [based  on  data  of  Sawaragi  and 
Deguchi  (1979)] 
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DIMENSIONLESS  LONGSHORE  TRANSPORT  RATE,  i' 


DIMENSIONLESS  DISTANCE  OFFSHORE,  x* 


Figure  (2.6).  Sediment  transport  and  velocity  profiles 

for  different  values  of  the  velocity  ratio 
u  . 

where  u  ,  is  the  oscillating  water 
w  mo 

velocity  magnitude  at  the  breaker  line  and 
w  is  the  sediment  fall  velocity.  [from 
Bailard,  1982] 
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DIMENSIONLESS  LONGSHORE  CURRENT.  V 


TABLE  I 


Location  of  Maximum  Longshore  Current  (v/v^ax) 
Longshore  Sand  Transport  (q/qnvax)  Distributions  From 

Different  Models 


Model  and 

Year 

x/x^^ 

Remarks 

Thornton 

(1973) 

— 

1.0 

Field  data  and 
theoretical  relation¬ 
ship  (2.11)  with 

B  =0.08 
s 

Komar 

(1977a) 

0.65 

0.80 

Theoretical  rela¬ 
tionships  (2.5), 

(2.12) 

Sawaragi 

and 

Deguchi 

(1979) 

0.40 

0.60 

Laboratory 

measurements 

Bailard 

(1982) 

0,60 

0.90 

_ 

Theoretical  rela¬ 
tionships 

_ 1 
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III.  DESCRIPTION  OF  LEADBETTER  BEACH  EXPERIMENT, 

SANTA  BARBARA,  CALIFORNIA 

A.  GENERAL 

Leadbetter  Beach,  Santa  Barbara,  California  was  chosen  as 
the  second  NSTS  experimental  site,  after  Torrey  Pines  Beach, 
California,  primarily  because  it  is  adjacent  to  a  total  sand 
trap  forming  a  spit  off  the  Santa  Barbara  breakwater,  and 
because  it  has  relatively  straight  and  parallel  nearshore  depth 
contours.  The  NSTS  Santa  Barbara  experiment  began  in  October, 
1979  with  a  14-month  sediment  trap  study  measuring  sand 
accumulation,  coupled  with  a  one-month  intensive  experiment 
during  27  January  through  25  February  1980. 

Exceptional  storms  were  encountered  towards  the  latter 
part  of  the  intensive  experiment.  During  the  period  16-20 
February  1980,  the  weather  conditions  deteriorated  greatly 
approaching  the  50-year  storm.  Significant  wave  heights  in 
excess  of  2  meters  and  longshore  current  speeds  of  150  cm/sec 
were  recorded.  These  conditions  coupled  with  spring  tides, 
high  winds  and  wave  set-up,  eroded  the  foreshore  and  back 
beach  areas  approximately  2  meters  in  elevation.  These 
phenomena  coupled  with  high  quality  data,  make  the  Santa 
Barbara  experiment  very  unique  and  significant  for  near-shore 
transport  studies  (Gable,  1981) .  This  study  emphasizes  the 
very  intensive  storm  period  of  17-18  February,  1980. 
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B.  COASTAL  GEOMORPHOLOGY 

The  shoreline  between  Point  Conception  and  Santa  Barbara 
has  the  unusual  (for  the  California  Coast)  east-west  orienta¬ 
tion  on  a  predominantly  north-south  coast.  The  coastal  areas 
of  Santa  Barbara  are  composed  of  marine  sedimentary  rock  and 
Monterey  Shale-bluffs  30  m  high  fronted  by  sand  and  cobble 
beaches  varying  from  0  to  30  m  in  width  (Corps  of  Engineers, 
1970) .  The  major  sources  of  sand  for  the  Santa  Barbara 
beaches  are  the  Santa  Maria  and  Santa  Ynez  rivers  (see  Fig. 
3.1) .  Some  sand  comes  from  the  erosion  of  the  coastal  cliffs 
(Trask,  1952) .  Bascom  (1951)  suggests  some  contribution  by 
wind  blown  sand,  but  neither  of  these  two  sources  are  rela¬ 
tively  importcint.  Textural  analysis  of  the  sands  indicates 
the  beaches  are  composed  of  fine  to  medium  grain  size  and  are 
relatively  well  sorted  with  a  mean  diameter  of  0.236  mm. 

C .  WEATHER 

The  climate  in  the  Santa  Barbara  area  is  classified  as  a 
Mediterranean  type  and  characterized  by  warm,  dry  summers  and 
mild,  wet  winters.  The  summers  are  mostly  dry  because  a 
semi-permanent  high  pressure  area  covers  the  Eastern  North 
Pacific  Ocean,  which  deflects  eastward  storms  to  the  north. 
During  winters,  this  Pacific  high  migrates  southward  and 
weakens,  allowing  occasional  frontal  systems  to  move  through 
Southern  California. 

The  prevailing  winds  for  the  West  Coast  are  from  the 
northwest  due  to  the  semi-permanent  Pacific  high  pressure  and 
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blow  generally  parallel  to  the  coastline.  However,  due  to 
the  orientation  of  the  coastline  at  Santa  Barbara  and  the 
protection  from  Point  Conception,  the  winds  are  generally 
light  and  variable  from  the  southerly  direction.  The  aver¬ 
age  wind  speed  is  less  than  3  m/sec  and  the  typical  pattern 
is  that  of  a  breeze. 

Most  extratropical  winter  storms  that  affect  the  Santa 
Barbara  area  move  southeastward  from  the  Northeast  Pacific; 
on  the  average,  1  or  2  such  storms  influence  the  Santa 
Barbara  channel  region  every  year. 

In  February  1980,  the  Eastern  Pacific  high  pressure  sys¬ 
tem  was  weaker  than  normal  and  displaced  southeast  of  its 
normal  position.  At  the  same  time,  the  low  pressure  center 
was  also  displaced  south  of  its  normal  position.  The  surface 
pressures  were  as  much  as  18  millibars  below  normal  east  of 
the  low  center.  The  Eastern  Pacific  high  pressure  system 
was  almost  non-existent,  which  shifted  the  storm  track  suffi¬ 
ciently  to  bring  storms  onshore  in  Central  and  Southern  Cali¬ 
fornia  (Miller,  1980)  . 

Two  series  of  frontal  systems  hit  the  study  site.  A 
strong  southeasterly  occurred  from  13  to  16  February  and  a 
very  intensive  southwesterly  from  17  through  21  February. 
Because  these  storms  were  southerly,  they  produced  large 
westerly  ocean  swells  and  waves  which  moved  into  the  Santa 
Barbara  channel.  The  relatively  warm  water  was  a  major  factor 
which  intensified  the  storm  and  brought  much  heavier  rain 
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than  normal.  Precipitation  was  13.2  cm  above  normal  at  this 
time  of  the  year.  All  these  factors,  coupled  with  extreme 
high  tides  and  wave  height,  caused  extreme  beach  erosion. 

D.  WAVES  AND  WIND  CLIMATE 

The  wave  climate  at  the  experimental  site  is  either  locally 
generated  inside  the  Channel  Islands  or  propagated  from  the 
open  ocean.  Normally,  the  predominant  waves  are  swell  from 
the  west-southwest  and  the  west  {240®-270°  true)  entering 
through  a  narrow  window  between  San  Miguel  Island  and  Point 
Conception  with  an  average  wave  height  of  91  cm  and  an  average 
wave  period  of  12  seconds  (Weigel,  1950).  The  waves  must 
refract  almost  90®  to  reach  Leadbetter  Beach  (Fig.  3.2).  As 
a  consequence,  the  swell  waves  reaching  Santa  Barbara  are 
greatly  reduced  in  height  and  arrive  at  large  angles  due  to 
the  severe  refraction  (O'Brien,  1950). 

As  a  result  of  the  origin  of  the  waves,  the  offshore 
islands,  and  the  refraction  of  the  waves,  the  wave  pattern 
in  Santa  Barbara  is  nearly  the  same  for  all  times  of  the  year, 
with  the  exception  of  the  waves  from  the  occasional  southeast 
storms.  The  swell  from  the  open  ocean,  refracting  and  break¬ 
ing  at  high  angle  of  incidence,  create  a  strong  easterly 
longshore  current.  The  occasional  southeasterly  storms  also 
create  reversals  where  the  longshore  current  is  to  the 
west. 

The  local  wind  waves  from  the  south  are  usually  insignifi¬ 
cant  with  only  choppy  seas  and  small  waves.  The  local  storms 
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Figure  (3.2).  Refraction  diagrams  for  Santa  Barbara  and 
vicinity  for  waves  of  10-second  period: 

(a)  waves  from  the  southwest;  (b)  waves  from 
the  south  and  southeast;  (c)  expanded  refrac¬ 
tion  drawing,  waves  from  the  west,  west- 
southwest,  and  southwest;  (d)  expanded  refrac¬ 
tion  drawing,  waves  from  the  southeast  (after 
O'Brien,  1950). 
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from  the  southeast  (125°-145‘’  true)  have  a  maximum  fetch  of 
145  Km  and  can  propagate  waves  toward  the  Santa  Barbara  Coast 
with  significant  wave  height  ranging  up  to  488  cm  (Weigel, 
1950) .  During  the  storm  period  (16-20  February,  1980)  ,  waves 
generated  thousands  of  ]cilometers  away  with  unlimited  fetch 
were  funneled  through  the  wave  window  of  (240°-270°  true) . 

The  wave  period  was  14-16  sec  with  significant  wave  height 
approaching  200  cm  from  the  southwest  and  west.  The  signifi¬ 
cant  wave  height  was  190  cm  from  the  southeast  on  February  16 
and  185  cm  from  the  west  on  February  20,  1980. 

E .  EXPERIMENT 

The  start  of  the  experiment  was  timed  such  that  a  strong 
negative  low  tide  occurred  during  the  daylight  hours  allowing 
nearshore  instruments  to  be  installed  on  a  dry  beach.  The 
tides  are  semidiurnal  with  a  mean  range  of  110  cm  referred 
to  mean  lower  low  water.  The  tides  ranged  from  -10  cm  to 
+200  cm  (MLLW)  during  the  intensive  experiment. 

A  plan  view  of  Leadbetter  Beach  (Fig.  3.3)  shows  the  loca¬ 
tion  of  the  instruments,  the  experiment  baseline  and  the  five 
range  lines.  It  is  noted  that  the  five  range  lines  are  not 
equally  spaced. 

Series  of  current  meters,  pressure  sensors  and  run-up 
meters  were  used  to  measure  the  wave  characteristics  during 
shoaling,  breaking,  and  run-up.  The  location  of  the  two  wave 
slope  arrays  is  shown  in  (Fig.  3.4). 
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Figure  (3.3).  Plan  view  of  the  study  site  s 
the  experimental  baseline  and 
five  ranges. 


1 .  Incident  Wave  Field  Measurements 

Two,  6  meter  square  arrays,  composed  of  four  Kulite 
solid  state  pressure  sensors  were  used  to  measure  an  un¬ 
biased  estimate  of  the  spectrum  of  the  radiation  stress, 

The  square  arrays  are  less  efficient  than  a  linear  array  in 
estimating  a  complete  directional  spectrum  but  the  narrow 
wave  window  provided  by  the  offshore  islands  of  Santa  Barbara 
Beach  produces  an  acceptable  wave  climate  estimate.  The 
location  and  orientation  of  the  two  slope  arrays  are  shown 
in  Fig.  (3.4).  The  water  depth  at  both  arrays  is  approximately 
9  m  and  the  sensors  were  mounted  120  cm  above  the  bottom. 

2 .  Profile  Measurements 

Profile  measurements  at  Santa  Barbara  consisted  of 
daily  wading  profiles  of  the  beach  face  and  shallow  nearshore, 
coupled  with  offshore  surveys  to  a  depth  of  4  m  using  a  trac¬ 
tor  and  boat  survey  for  pre-  and  post-experiment  measurements 
to  12  m  depth.  Five  range  lines,  extenuing  essentially  per¬ 
pendicular  to  the  beach  and  covering  200  m  alongshore,  were 
established  for  the  daily  wading  profile  (Fig.  3.3).  The  off¬ 
shore  distances  along  range  lines  are  referenced  to  the  distance 
from  the  experimental  baseline  (y-axis) . 

The  wading  profiles  were  done  using  conventional  rod 
and  level  survey  techniques  on  a  daily  basis.  Range  lines 
were  defined  by  two  survey  points  marked  by  flags,  which  a 
rodman  aligned  before  each  survey  reading.  Rod  stations  were 
measured  using  a  specially  constructed  plastic-coated  steel 
survey  line  marked  at  5  m  intervals.  An  engineer's  level  and 
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rod  provided  vertical  control.  The  profile  survey  was  ter¬ 
minated  when  the  water  became  too  deep  for  the  rodman,  or 
the  breaking  waves  made  it  impossible  to  plumb  the  rod.  These 
surveys  were  completed  at  the  lowest  possible  tides  occurring 
during  daylight  hours.  Each  svirvey  was  closed  to  the  original 
benchmark  or.  land  to  verify  survey  accuracy.  Data  were 
entered  on  site  by  hand  into  a  Tektronix  computer  terminal 
within  a  short  time  after  each  survey.  Profiles  were  plotted 
immediately  to  visually  search  for  errors  which  were  due  to 
either  data  entry  or  a  bad  survey  point  which  could  be  cor¬ 
rected  by  redoing  the  field  survey  immediately. 

Generally,  small  scale  changes  occurred  in  the  profile 
during  most  of  the  intensive  month-long  experiment  (Fig.  3.5), 
except  during  the  storm  where  the  beach  was  cut  vertically 
2.5  m  and  horizontally  about  60  m  within  almost  a  four  day  period. 
Although  it  was  impossible  to  maintain  the  nearshore  instru¬ 
mentation  during  the  storm,  the  deep  water  instruments  con¬ 
tinued  to  work,  providing  an  extremely  valuable  data  set  on 
beach  erosion. 
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IV.  FIELD  DATA  ANALYSIS 


A.  BEACH  PROFILES  AND  CONSERVATION  OF  SAND  TRANSPORT  EQUATION 
Five  beach  normal  ranges  at  Leadbetter  Beach  were  surveyed 
on  a  daily  basis  starting  from  the  established  baseline  on 
the  beach  to  a  distance  about  100  m  offshore.  Standard  rod 
and  level  technique  described  in  (III.E.b)  were  used.  Data 


from  the  boat-fathometer  survey  conducted  on  25  February,  1980 
for  the  five  ranges  were  used  to  extend  each  profile  out  to 
the  shallow  water  limit  at  x  =  180  m  where  ^  -^) . 

Linear  interpolation  was  applied  twice,  first  to  calculate 
the  missing  data  points,  and  second  to  define  five  equally 
spaced  ranges  (50  m  apart) .  The  five  ranges  were  averaged 
for  both  17th  and  ISth  of  February  to  produce  a  daily  mean 
profile.  Then,  a  "mean  beach"  prof ile  was  defined  by  averaging 
these  two  daily  profiles.  Fig.  (4.1)  shows  the  mean  beach 
profile,  relative  to  mean  sea  level  (MSL)  and  indicates  no 
well  developed  bar  structures  or  trough  along  the  profile. 


Based  on  the  difference  in  elevations  between  the  17th 


and  18th  February  surveys'  at  each  location,  the  net  longshore 
transport  distribution  is  calculated  using  the  conservation 
analysis  of  sand  transport  described  later  in  this  section. 

A  "mean  net  longshore  transport",  distribution  was  ob¬ 

tained  by  averaging  the  net  transport  distributions  for  the 
five  ranges.  The  "mean  net  sand  transport"  distribution  and 
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the  standard  deviation  calculated  at  each  point  are  shown  in 

Fig.  4.2  indicating  a  small  variance  with  acceptable  scatter 

around  the  q  distribution. 

^mean 

The  q  ^  „  distribution  goes  to  zero  after  a  distance  of 
^mean  ^ 

100  m  offshore  as  shown  in  Fig.  4.2  because  the  wading  pro¬ 
files  on  the  17th  and  18th  are  extended  after  this  distance 
by  the  same  fathometer  data  taken  on  25  February. 

Profile  changes  between  sequential  measurements  can  be 
approximated  by  a  parallelepiped  with  sides  Ax,  Ay  and  the  time 
varying  depth  (Ah) .  Considering  the  equilibrium  of  the 
parallelepiped  element  (Fig.  4.3),  the  net  sand  trnasport  rate 

in  the  x-direction  is  q  4.  Ay  =  — ^AxAy  and  the  net  sand 

^x  net  3x  ~ 

transport  rate  in  the  y-direction  is  q  Ax  =  -:^AyAx. 

The  total  net  sand  transport  rate  (q^  Ay  +  q^  Ax)  is 
balanced  by  the  rate  of  the  net  volume  changes  over  the  ele¬ 
ment  giving 


'^xnet'^y  Snet^^ 


(1-a)  ^  AxAy 


(4.1) 


where  a  is  the  sand  porosity  and  is  taken  as  0.40.  Equation 
(4,1)  can  be  written  as 


3x  3y 


, ,  .  3h 

3t 


Since  a  strong  longshore  current  was  observed  during  17 


and  18  February,  (see  Fig.  (6.3)),  a  significant  longshore 
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sand  transport  was  expected.  Therefore,  it  is  hypothesized 
that  the  longshore  sand  transport  is  much  larger  than  the  on- 


offshore  sand  transport,  i.e.,  >> 


The  profile  changes 
in  Fig.  3.6  do  not  show  any  well-defined  bar  or  trough 


structures  which  suggests  the  on-offshore  sand  transport  is 


negligible  (i.e.. 


9x 


=  0).  Thus,  equation  (4.2)  becomes 


3q 

Ty 


=  (1-a) 


3h 

at 


(4.3) 


Equation  (4.3)  states  simply  that  the  changes  in  the  beach 
profile  with  respect  to  time,  are  proportional  to  the  net 
longshore  sand  transport.  The  transport  rate/unit  width  of 
on-offshore  distance,  q^,  can  be  calculated  since  water  depth 
(h)  is  measured  as  a  function  of  time  (t) .  The  boundary  con¬ 
ditions  are  needed  to  solve  (4.3).  For  the  offshore  boundary 
condition  it  is  assumed 


q  =  0  at  h  =  h 

^y  o 


where  h  is  the  water  depth  at  the  shallow  water  limit, 
o 

h^  =  and  L  is  the  wavelength.  Given  the  dominant  wave 

period  of  14  sec  (f  =  0.0703  Hz) ,  h  was  calculated  to  be 
5.0  m  at  a  distance  approximately  180  m  from  the  baseline. 

The  nearshore  boundary  condition  requires  the  sand  transport 
to  go  to  zero  at  the  shoreline.  It  appears  from  Fig.  4.2 
that  this  condition  is  satisfied  from  the  data  at  the  baseline 
at  X  =  0.  The  net  sediment  transport  rate  is  calculated  by 
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integrating  the  net  cross-shore  distribution  which  is  assumed 
to  be  similar  to  the  longshore  sand  trauisport  distribution. 

For  discussion  purposes,  two  regions  within  the  surf  zone 
are  defined.  A  swash  region  is  defined  from  the  baseline  to 
the  intersection  of  HWL  with  the  profile,  i.e.,  (from,  x  =  0 
to  X  -  30  m) ;  a  seaward  region  is  defined  from  x  -  30  m 
to  X  -  180  m.  The  swash  region  indicates  significant 
erosion  on  the  beach  face  (Fig.  4.2),  which  will  not  be  considered 
in  this  study  because  the  model  does  not  include  sediment 
transport  on  the  beach  face.  The  seaward  region  shows  the 
eroded  profile  with  a  small  accreted  area  around  x  =  90  m, 
which  appears  with  a  negative  sign. 

B.  RADIATION  STRESS,  WAVE  HEIGHT  AND  DIRECTION  AT  OUTER 

BOUNDARY 

Measurements  from  a  610  cm  square,  four  pressure  sensor 
array  in  8-m  water  depth  are  used  to  produce  an  unbiased 
estimate  of  the  spectrum  of  radiation  stress.  The  pressure 
arrays  worked  throughout  the  experiment  including  the  storm 
period.  The  array  measurements  are  also  used  to  calculate 
the  height  cind  the  directional  characteristics  of  the  waves. 

The  radiation  stress  (S  )  is  given  by  (Longuet-Higgins , 

1962) 


n 

S  =  /  pu^(z,t)v  (z,t)  dz  (4.4) 

xy  w 


where  u  and  v  are  the  crossshore  and  alongshore  components 
w  w 
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of  wave  velocity.  The  velocities  are  calculated  from  the 
pressure  gradients  measured  at  the  array  applying  linear  wave 
theory  spectral  treinsf ormations.  For  a  particular  component 


Uj (f ,z,t) 


1 fi  cosh  k (h  +  2) 1  9p  ■  _  t  5 

0)  p  cosh  k  (h  +  z  T  3  X .  '  ^  ' 

m  j 


(4.5) 


where  the  z  is  the  measured  depth  of  the  pressure  sensors, 
m 

oj  is  the  radial  frequency  and  j  refers  to  the  two  horizontal 
components.  The  term  between  the  brackets  is  the  spectral 
transfer  function  H(f)  relating  the  measured  pressure  gradients 
to  inferred  velocity  components.  Substituting  equation  (4.5) 
into  (4.4),  using  the  measured  pressure  head  ^  and  integrating 


-]  C, 


S  (f)  =  [ - 5 - j  ^ 

^  2  k(cosh^k(h  +  z  )) 

m  o X  <3y 


(f) 


(4.6) 


where  co-spectrum  between  the  measured 

lf'3y 

pressure  gradient  components. 

Wave  heights  were  determined  from  the  surface  elevation 
records  using  the  "zero-up-crossing  method"  in  which  the  wave 
height  is  defined  as  the  difference  between  the  maximum  and 
minimum  occurring  between  two  consecutive  zero-up-crossings. 
Surface  elevation  records  were  first  linearly  detrended  to 
exclude  the  effects  of  the  rising  and  falling  of  the  tides 
and  then  high-pass  filtered  with  a  cut-off  frequency  of  0.05 
Hz  to  exclude  surf  beat,  followed  by  a  low-pass  filter  with 
a  high  frequency  cut-off  (0.3  Hz). 
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The  filtering  was  accomplished  by  Fourier  transforming  the 
signals,  zeroing  the  Fourier  cunplitude  coefficients  in 
the  filtered  frequencies,  and  inverse  transforming  the 
complex  spectrum  to  obtain  the  filtered  time  series.  The 
entire  68-minute-record  was  transformed  at  one  time  to  mini¬ 
mize  the  end  effects  which  result  in  spectral  leakage,  and 
to  obtain  maximum  resolution  giving  very  sharp  roll-off  at 
the  filter  cut-offs. 

Due  to  the  large  number  of  current  meters ,  current  data 
are  utilized  applying  linear  theory  to  infer  wave  heights 
from  the  velocities  (Guza  and  Thornton,  1980) .  Surface  ele¬ 
vation  time  series  were  obtained  from  current  measurements 
using  the  complex  Fourier  spectra  of  the  horizontal  velocity 
component  U(f),  V(f).  The  complex  surface  elevation  spectrum, 
X(f),  is  calculated  applying  the  linear  wave  transfer  function, 
H(f) 


X(f)  =  H(f)  •  V(f)  (4.7) 

Then,  the  complex  surface  elevation  spectrum  is  inverse 
transformed  to  obtain  the  svurface  elevation  time  series  from 
which  the  wave  height  distribution  is  calculated.  Surface 
elevations  were  also  inferred  from  pressure  signals  in  a  simi¬ 
lar  manner  by  transforming  the  pressure  records  using  linear 
theory. 

A  mean  incident  wave  angle  (a)  and  a  representative  fre¬ 
quency  (fp)  are  defined  for  input  into  the  models  such  that 
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using  a  and  f  to  calculate  the  total  radiation  stress  S  , 
P  xy' 

the  same  value  is  obtained  as  that  measured  by  the  pressure 
sensor  array. 

The  radiation  stress  (S^^)  in  equation  (4.4)  can  be  ex¬ 
pressed  applying  linear  wave  theory  for  long-crested  waves 


S^y(f)  =  E(f)  n(f)  sin  a  cos  a  (4.8) 


The  total  radiation  stress  (S  i  is  calculated  using 

xy  ’ 


S  =  /  S  (f)  df 

xy  •’  xy 


(4.9) 


T  T 

Then,  the  calculated  is  equated  to  its  equivalent 

in  (4.8)  for  monochromatic  (f  =  f  )  waves  approaching  at  a 

It 

mean  angle  a 


S  n(f  )  sin  a(f  )  cos  a(f  )  (4.10) 

xy  P  P  P 


T 

where  E  is  the  total  energy.  Applying  linear  theory  rela¬ 
tionships  a  mean  incident  wave  angle  a  is  defined  as: 


a(fp) 


1 

I 


sin 


-1 


(4.11) 


where  the  representative  frequency  f  is  the  peak  frequency 

P 

of  the  energy  spectrum.  Since  the  radiation  stress  is  calcu¬ 
lated  using  the  pressure  gradients,  the  kinetic  energy  calculated 
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using  the  pressure  gradient  is  introduced  for  consistency 
of  the  analysis.  It  is  assumed  the  kinetic  energy  is 
half  the  total  energy.  Then  equation  (4.11)  yields 


The  kinetic  energy  is  calculated  using  the  pressure  gradient 
component  spectrum 


Ej^(f)  =  !H{f)l^[Gg  (f)  +G5  (f)]  (4.13) 

3x  3y 

I  I  2 

where  |H{f) |  is  the  linear  wave  transformation  function 
given  by 

!H(f)i'^  =  — sinh(2kh) (4.14) 

4ka)'^  cosh"^k(h  +  z  ) 
m 

The  mean  incident  wave  angle  (a)  at  the  pressure  sensor's 
location  was  then  refracted  shoreward  to  the  5  m  contour  (at 
the  shallow  water  limit) .  A  refraction  calculation  was  per¬ 
formed  over  a  bathymetric  chart  of  the  experimental  site  that 
shows  the  beach  has  straight  and  parallel  contours  shoreward 
of  5  m  depth  (Fig.  4.4).  Applying  linear  refraction  to  the 
measured  mean  angle  (a  =  12.4°)  at  the  pressure  sensor  location 
and  peak  frequency  (fp  =  0.07  Hz),  waves  reach  the  5  m  contour 
with  a  refracted  angle  =  9.0°  where  the  calculated  refrac¬ 
tion  coefficient  equals  0.985. 
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The  is  calculated  using  linear  refraction  theory: 


H 


rms 


H 


=  ( 


rms  (8m) 


<-r> 


=  (45=) 


1/2 


(Kr) 


(4.15) 


giving  a  refracted  equivalent  to  .85  m  at  water  depth 

h^  =  5  m,  which  is  used  as  the  initial  wave  height  for  the 

model.  Refracted  wave  rays  were  extended  to  get  the  angle 

of  approach  at  the  breaker  line  calculated  to  give  =  6.0°. 

The  breaker  line  was  defined  at  maximum  H  .  a.  was  used 

rms  b 

to  calculate  the  longshore  sand  transport  using  the  wave 
power  expression. 

Table  II  gives  various  calculated  parameters  of  a, 

f_,  H  „  at  different  times  before,  after  and  during  the 
p  rms  ^ 

17th  and  18th  of  February.  Orientation  of  the  coordinates 
leads  to  a  negative  sign  for  both  and  a,  indicating  the 

waves  are  coming  from  the  west. 
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TABLE  II 


Calcxilated  Parameters  at  West  Pressure  Sensor 
(Diiring  17-18  February) 


1  Date 

time 

depth 

(m) 

tide 

(cm) 

(Erg/cm  ) 

(Hz) 

rms 

(cm) 

o 

00 

(ys 

rH 

0730 

9.2 

+71 

-0.43x10^ 

-5.8 

0.094 

i 

65.6 

1 

5 

1800 

7.8 

-71 

-1.76x10^ 

-12.4 

0.070 

00 

00 

00 

rH 

2230 

9.4 

+  86 

-2.16x10^ 

■■ 

-13.0 

0.062 

96.2 

O 

CO 

cr» 

rH 

0920 

9.5 

+96 

-2.61x10^ 

-15.9 

0.070 

96.2 

F 

1 

2200 

_ 

9.2 

+  74 

5 

-1.25x10^ 

-15.0 

—  . 

0.094 

o 

o 
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V.  LONGSHORE  SEDIMENT  TRANSPORT  MODEL 


A  littoral  sand  transport  model  is  developed  based  on 
the  energetics  approach  expressed  by  Thornton  (1973)  in  which 
the  energy  transformation  is  described  by  the  random  wave 
model  of  Thornton  and  Guza  (1983a)  coupled  with  the  longshore 
current  model  of  Thornton  and  Guza  (1983b) .  The  input  to 
the  predictive  model  are  the  offshore  rms  wave  height,  mean 
direction  and  mean  period.  Both  an  analytical  model  for  plane 
sloping  beaches  and  a  nimerical  model  for  arbitrarily  bottom 
profile  are  presented. 

A.  DEVELOPMENT  OF  SAND  TRANSPORT  MODEL 

Assuming  stationary  wave  conditions  with  a  beach  of 

straight  and  parallel  contours,  the  gradient  of  the  wave 

dECg 

power  in  the  on-offshore  direction  — is  balanced  by  the 
ensemble  averaged  dissipation  due  to  both  wave  breaJcing 
<£j^>  and  bottom  friction  <£f>  (Thornton  and  Guza,  1983)  such 
that 


dECg^ 

dx 


<e  >  +  <t 

D 


(5.1) 


The  wave  power  is  expressed  using  the  linear  wave  theory 
relationships : 


^  ^  P  5  “Ls  =  P(H)  dH  (5.2) 

o 
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2  k  h 


]  cos  ct 


(5.3) 


= 


C 

2 


[1  + 


sinh  ikh 


where  is  the  root  mean  square  wave  height,  a  is  the 

mean  wave  direction  and  k  is  the  wavenirrher  associated  with 
peak  frequency  f^. 

Dissipation  due  to  wave  breaking  can  be  modeled 

after  a  periodic  bore,  applying  conservation  of  mass  and  momen¬ 
tum  at  regions  of  uniform  flow  up-and-downstream  of  the  bore 
as  shown  in  Fig.  (5.1).  The  average  rate  of  energy  dissipa¬ 
tion  per  unit  area  for  single  breaking  wave  is  calculated 
(Stoker,  1957) : 

2  (h_  -  h,  )  ,  4 

We  =  4^5  (5.4) 

12  h 

where  the  wave  height  H  is  measured  as  the  maximum  to  mini¬ 
mum  of  the  bore,  Q  is  the  volume  discharge  of  flow  per  unit 
area  across  the  bore  described  by  simple  linear  bore  theory 
after  Hwang  and  Divoky  (1972) ,  and  B  is  a  breaker  coefficient 
of  0(1).  The  coefficient  B  is  an  lanspecified  parameter  in 
the  model  ,  which  is  determined  from  the  data. 

The  Rayleigh  distribution  is  assumed  to  describe  the 
waves  everywhere  with  the  implied  assumptions  that  the  waves 
are  very  narrow  band  such  that  they  can  be  described  by  a 
single  frequency  f  and  mean  direction  a.  Thornton  and  Guza 

Sr 

(1983)  predict  the  fraction  v/aves  which  are  breaking  at 

each  location.  Using  this  information,  the  average  rate  of  energy 
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Figure  (5.1).  Periodic  bore  used  to  describe 
spilling  breakers. 
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(5.5) 


dissipation  due  to  breaking  <ejj>  is  given  by 

S/tt 


B^f 


<eb>  = 


16 


p  g 


_ E 

4..  5  rms  ' 


Y  h 


where  y  is  a  coefficient  relating  H  to  h  (taken  to  be  0.44 

rms 

by  Thornton  and  Guza,  1983) .  The  average  frictional  dissi¬ 
pation  <££>  at  the  bottom  boundary  layer  for  the  ensemble 
(all  waves)  calculated  using  Rayleigh  probability  distribution 
is : 


P 

16/? 


2tt  f  H 
j.  p  rms. 

sinh  kh  ^ 


(5.6) 


where  is  the  friction  coefficient.  Thornton  and  Guza 
(1983a)  point  out  that  <£^>  is  relatively  more  important  for 
gently  sloping  beaches  and  for  low  frequency  waves,  i.e.,  the 
shallower  the  depth  the  greater  the  frictional  dissipation. 
They  neglect  the  frictional  dissipation  as  it  has  a  small  con¬ 
tribution  compared  with  as  will  be  done  here  also. 

Substituting  equations  (5.2),  (5.3),  (5.5),  and  (5.6) 
into  equation  (5.1) ,  yields: 


d  1 
dx  S’ 


p  g  H 


2 

rms 


C  cos  a 

g 


3^  p  g  7 

16  4.5  "rms 

Y  h 


(5.7) 


^rms  evaluated  analytically,  assuming  small  inci¬ 

dent  wave  angle  (cos  a  -  1)  and  limiting  the  analysis  to 
shallow  water  (^  ^  ^)  ,  which  gives 
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(5.8) 


F 


H 


rms 


=  a 


1/5 


h9/10ji.h23/4^. 


1 

T377 


)  ] 


-1/5 


where 


and 


y  = 

-'o  o  o 


The  value  of  tan  8  is  the  beach  slope  and  the  siibscript  "o" 
refers  to  the  input  conditions  at  the  outer  shallow  water 
depth  limit.  When  the  water  depth  gets  very  shallow,  (5.8) 
become  s 


H 


rms 


=  a 


1/5  j^9/10^. 


as 


(5.9) 


which  indicates  that  the  wave  height  in  the  inner  surf  zone 
is  related  to  the  depth  and  independent  of  the  initial  condi¬ 
tions  in  deeper  water. 

With  the  assumptions  of  straight  and  parallel  depth  contours 
and  stationary  wave  conditions,  the  longshore  current  dis¬ 
tribution  across  the  surf  zone  can  be  estimated  using  the 
longshore  momentum  flux  balance  (Thornton  and  Guza,  1983) 


3S 


3x 


-  T, 


(5.10) 
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is  the  radiation 


where  x,  is  the  bottom  shear  stress  and  S 

b  xy 


stress  defined  by  Longuet-Higgins  and  Stewart  (1962)  (for 
linear  theory)  as 


xy 


or-  «  Sin  a 

E  Cg  cos  a  — g — 


From  linear  wave  refraction 


sin  CL 


sin  a 


=  constant 


Thus , 


as. 


JiZ  = 


sin  a 


ax 


“  ^  ^g 


(5.11) 


The  term  in  brackets  in  the  R.H.S.  is  defined  in  (5.1) 
as  <Cj3>.  Therefore  (5.11)  reduces  to; 


sin  a 

o 

C 

o 


"b 


(5.12) 


The  average  alongshore  bottom  shear  stress  Xj^  is  defined  by 
Longuet-Higgins  (1970) ,  assuming  the  quadratic  law 


Tfa  =  0  Cg  ju^l (V  +  v^): 


(5.13) 


where  u^  is  the  total  instantaneous  velocity  due  to  waves 
and  current,  V  is  the  longshore  current,  v^  is  the  longshore 
component  of  the  oscillating  wave  velocity  and  c^  is  the  bed 
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shear  stress  coefficient.  The  bar  over  equation  (5.13)  indi¬ 
cates  averaging  in  time  over  one  wave  period.  Assuming  a 
small  angle  of  wave  incidence  (a  <  15®) ,  and  a  weak  current 
case  following  Liu  Ph.  L-F  and  R.  Dalrymple  (1978) ,  equation 
(5.13)  simplifies  to 


where  u^  is  the  mean  cross-shore  wave  velocity. 

I u^l  is  specified  using  the  shallow  water  wave  theory 
relationships  and  Rayleigh  wave  height  distribution  (Thornton 
cind  Guza,  198  3)  such  that: 


w' 


T  "rms 


(5.15) 


Substituting  (5.11),  (5.12),  (5.14),  (5.15)  into  (5.10) 
and  solving  for  the  mean  longshore  current  velocity,  yields: 


V 


sin  a. 


P 


u 


w ' 


3  f  g^^^  sin  a 
3  p  ®  o  rms 

T  „  4  c  .9/2 

Cj  Y  oh' 


(5.16) 


The  longshore  sand  transport  distribution  is  obtained 

by  substituting  the  expressions  for  the  change  in  energy  flux 
3EC 

(— r— ^) ,  equation  (5.1)  and  the  longshore  current  equation 

d  X 

(5.16)  into  equation  (2.11)  to  give 
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q 


s 


(5.17) 


Substituting  (5.7),  (5.15)  and  (5.16)  in  equation  (5.17) 

yields: 


B 


g(l  -^) 


3  9  19 

sin  a  B  1/2 

o  p  rms  1  ' 

n  ^  i 


(5.18) 


Equation  (5.18)  has  the  observed  characteristics  that  outside 
the  surf  zone  as  h  reaches  h^,  qg(x)  gets  small  and  when  h 
goes  to  zero  at  the  shoreline,  qg(x)  approaches  zero.  A 
maximum  transport  occurs  between  these  two  boundaries. 

For  a  planar  sloping  beach,  the  analytical  solution  for 
H___  (5,8),  which  is  limited  to  shallow  water,  is  substituted 

X^IuS 

into  (5.18)  to  obtain  the  cross-shore  transport  distribution. 

For  arbitrary  bottom  profiles  (but  straight  and  parallel  depth 

contours) ,  a  simple  forward  differencing  scheme  is  used 

to  numerically  integrate  equation  (5.1)  for  H  .  The 

rms 

calculated  H  value  is  substituted  into  (5.18).  The  energy 
rms 

flux  balance  is  also  numerically  integrated  from  initial 
offshore  distance  (at  depth  h^)  to  the  shoreline  with  a 
spatial  interval  Ax 


E  Cg 


x'  i+1 


=  E  Cg 


X'  i 


-  Ax'|<ej^> 


i  * 


1  — "  1,2, ...,n 


(5.19) 
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Equation  (5.19)  starts  at  the  initial  offshore  distance 


(i  =  1)  where  H  ,  h_,  f_,  a  are  given.  The  model  calculates 
o  o  p  o 

the  forwaurd  step  by  subtracting  the  total  dissipation  from 
the  calculated  energy  flux  evaluated  from  the  previous  step. 
The  model  predicts  each  step,  using  linear  theory  rela¬ 

tionships  (5.2)  and  (5.3)  until  i  approaches  n  at  the  shore¬ 
line.  Then,  the  predicted  H  values  are  used  to  calculate 

•L  XllS 

the  longshore  current  using  (5.16).  The  sand  transport  dis¬ 
tribution  is  evaluated  using  (5.18). 

Integrating  q  in  (5.18)  over  limits  of  offshore  distances 
s 

for  the  study  zone,  the  volume  transport  rate  (Q  )  is  calcu¬ 
lated  as 


o  _ 

=  /  q^  dx  (5.20) 


where  x  is  the  offshore  distance  at  which  the  initial  condi- 
o 

tions  of  waveheight  (H^)  and  water  depth  (h^)  are  given. 

The  coefficient  B  in  (5.18)  is  used  to  calibrate  the  model. 

s 
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I.  MODEL  RESULTS  AND  COMPARISON  WITH  PROTOTYPE  DATA 


A.  MODEL  RESULTS  FOR  H _  AND  V 

rms 

Field  data  for  the  Santa  Barbara  experiment  are  utilized 

to  calculate  the  parameters  a^,  and  which  in 

turn  are  used  as  input  for  both  the  expressions  (5.8)  and 

(5.19).  Model  coefficients  y,  B  and  c,  were  determined  from 

H  ^ 

the  data.  The  ratio,  y  =  — ^  was  found  to  be  0.44  from  the 
slope  of  the  h  vs,  H  curve.  Fig.  (6.1) ,  in  the  inner  surf 
zone.  This  agrees  with  the  value  reported  earlier  by  Thornton 
and  Guza  (1983)  for  the  Torrey  Pines  Beach  experiment. 

^rms  inferred  from  the  pressure  and  velocity 

measurements  at  time  1800  on  the  17th  of  February  are  compared 
with  the  wave  transformation  model  for  various  breaker  coeffi¬ 
cients,  B,  Using  an  iterative  scheme,  based  on  the  least 
squares  error  criterion,  B  equal  to  1.1  gave  the  best  fit  to 

H  data  (Fig.  6.2a),  To  check  on  the  generality  of  the 
rms 

model,  the  waveheight  transformation  model  was  run  for  5  Feb¬ 
ruary  using  B  =  1.1,  It  also  gave  the  best  model  fit  to 
the  data  (Fig.  6.2b). 

An  optimal  bed  shear  stress  coefficient,  c^,  was  calculated 
comparing  the  longshore  current  model  with  the  measured  long¬ 
shore  current  values,  V.  Since  there  are  only  three  longshore 
current  measurements  available  on  17  February  due  to  the  des¬ 
truction  of  instruments  by  the  storm,  the  data  from  5  Feb¬ 
ruary  were  again  used.  With  B  =  1.1,  c^  was  varied  to 
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Figure  (6.1) .  Root  mean  square  waveheight  vs.  water 
depth  for  5  February,  1980,  giving 
Y  =  0.44  in  the  inner  surf  zone. 
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Model  H  as  a  fvinction  of  offshore  distance 
comparea  with  measurements  indicated  by  x.  T 
model  is  applied  to  data;  (a)  on  17  February; 
(b)  on  5  February,  to  obtain  B  which  was 
found  to  be  1.1  in  both  days. 


0.005 


obtain  the  best  fit  to  the  measured  V.  A  value  of  = 
gave  the  best  fit  (Fig.  6.3b). 

B.  MODEL  RESULTS  FOR  CROSS-SHORE  SAND  TRANSPORT  DISTRIBUTION 

The  empirical  coefficient  B^  in  equation  (5.18)  is  used  to 

calibrate  the  model  and  will  be  explained  later.  In  all  the 

calculations  and  figures,  the  offshore  distance  (x)  is 

measured  starting  from  the  experimental  baseline  at  x  =  0. 

The  average  net  transport  between  profiles  on  17  and  18 

February  are  to  be  determined.  A  weighted  time  difference 

of  25  hours  was  found  between  the  two  mean  profile  surveys. 

The  model  predictions  of  H  ,  V  and  q  at  each  x-location 
^  rms  ^s 

depend  mainly  on  the  water  depth  (h) .  Therefore,  two  cases 
are  considered  to  test  the  sensitivity  of  the  models  to  the 
variations  in  the  water  depth,  in  the  first  case,  the  water 
depth  does  not  vary  with  time.  In  the  second  case,  the  water 
depth  varies  by  hourly  time  increments  according  to  tidal 
variation. 

1 .  Model  Results  for  Constant  Water  Depth 

In  this  case  study,  the  water  depths  are  held  constant 

at  relative  MSL.  The  numerical  model  utilizes  the  actual 

bottom  and  predicts  H  ,  V  and  q  .  The  model  results  are 

rms  ^s 

svimmarized  in  Fig.  6.4  indicating  the  capability  of  the 
model  to  predict  the  increase  in  waveheight  due  to  shoaling 
and  the  decrease  due  to  breaking  processes.  The  longshore 
sand  transport  and  the  longshore  ciorrent  have  a  similar  dis¬ 
tribution,  which  implies  that  the  sand  transport  may  primarily 
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Figure  (6.3).  Longshore  current  (68-minute  average)  as  a 
function  of  offshore  distance  compared  with 


measurements  indicated  by  x.  The  model  is 
applied  to  data  on;  (a)  17  Feb.,  showing  only 
a  few  instruments  working;  (b)  5  Feb.,  for  which 
c-  =  0.005  was  found  to  give  the  best  fit. 


be  governed  by  the  longshore  current  (Fig.  6.4).  The  influ¬ 
ence  of  rapid  depth  changes  is  reflected  in  the  small  pertur¬ 
bation  of  the  V  and  curves  at  x  =  65  m,  which  is  associated 
with  a  local  flatness  of  the  beach  slope,  as  can  be  seen  in 
the  "mean  beach  profile,"  Fig.  4.1.  The  maximum  waveheight 
is  used  to  define  a  "mean  breaker  line"  location,  x^^,  for 
discussion  purposes  and  normalization  of  the  results.  The 
analytical  model  is  applied  to  a  plane  sloping  beach.  A  slope 
of  0.037  was  found  as  the  best  fit  to  the  "mean  beach"  pro¬ 
file  using  Fig.  4.1.  The  normalized  distributions  resulting 
from  the  analytical  model  are  summarized  in  Fig.  6.5,  giving 
the  maximum  sand  transport  at  x/x^  -  0.80,  which  is  similar 
to  the  results  of  Komar  (1977a) . 

Comparison  between  the  sand  transport  distributions  using 
the  numerical  and  analytical  models  is  given  in  Fig.  6.6. 

The  sand  transport  peak  of  the  numerical  model,  which 
uses  the  actual  bathymetry,  is  shifted  offshore  compared 
with  the  analytical  distribution  using  a  planar  sloping  beach. 
It  may  be  seen  in  equation  (5.18)  that  q^  is  proportional 

to  (-^)  ,  which  means  q  increases  as  h  decreases  for  the 

h '  ® 

same  location.  The  actual  bottom  profile  (Fig.  (4.1))  used 
in  the  numerical  model  is  deeper  near  the  breaker  line  com¬ 
pared  with  the  plane  sloping  beach  used  in  the  analytical 
model;  therefore,  the  calculated  maximum  q^  in  the  numerical 
model  is  shifted  offshore  relative  to  the  analytical  model. 
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Figure  (6.4).  Model  predictions  of  H  (••••)»  V  ( - ) 

and  q  ( - )  for  actui?®bathyinetry  holding 

depth^constant. 
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Figure  (6.5).  N^ormalized  distribution  for  H  ,  V  and 
q  calculated  by  the  analyticiT^model 

(constant  depth) ;  H  {••••),  V  ( - ) 

and  qg  ( - )  . 
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Figure  (6.6).  Average  sand  transport  distribution 

calculated  for  constant  depth  case;  for 

analytical  model  ( - )  and  numerical 

model  ( - )  . 
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2 .  Model  Results  for  Time-Varying  Water  Depth 

The  models  predicting  H  ,  V  and  q  were  run  for 
the  25  hours  using  a  one  hour  time  step  and  correcting  the 
depths  for  the  tidal  elevations  at  each  time  step.  The 
resulting  25  sand  trainsport  distributions  were  averaged  over 
the  25  hours  to  give  a  predicted  mean  longshore  sand  trans¬ 
port  distribution  between  profile  surveys  on  17  and  18 
February. 

Mean  sand  transport  distributions  including  tidal 
variation  of  both  the  analytical  and  numerical  models  are 
given  in  Fig.  6.7  which  shows  a  good  agreement  with  a  peak 
shift  towards  offshore  for  the  nvimerical  model  following  the 
previous  discussion.  When  considering  the  tides,  the  averaging 
causes  the  distribution  to  spread  onshore  to  x  «  30  m  with 
less  q„  „  than  the  first  case  for  the  same  coefficients. 

In  both  cases  the  maximum  longshore  sand  transport  occurs  at 
almost  the  same  location  as  the  longshore  current  which  indi¬ 
cates  again  that  the  sand  transport  is  primarily  driven  by 
the  longshore  component  in  both  models.  The  sediment  trans¬ 
port  distribution  predicted  by  the  numerical  model  is  not 
smooth,  unlike  the  analytical  curve.  This  is  because  of  the 
use  of  the  actual  depth  profile. 

C.  MODEL  RESULTS  AND  COMPARISON  WITH  PROTOTYPE  FIELD  DATA 
The  net  longshore  transport  rate  distribution  for  Santa 
Barbara  Beach  is  calculated  using  the  profile  changes  during 
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Figure  (6.7).  Average  sand  transport  distribution  over 
the  daily  tidal  cycle  calculated  using 
variable  depth;  for  analytical  model 
( - )  and  nvanerical  model  ( - )  . 
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17  and  18  February.  Most  of  the  estimated  sediment  trans¬ 
port  occurs  at  a  distance  20  m  from  the  baseline  indicating 
a  significant  beach  face  erosion,  which  is  beyond  the  scope 
of  the  present  study. 

Since  the  volume  sand  transport  rate  calculated  from  the 
data  is  the  net  volume  chcinge  rate  in  the  beach  profile,  and 
does  not  represent  the  total  longshore  sand  transport  as  the 
model  predicts,  care  must  be  taken  in  the  interpretation  of 
the  comparisons. 

A  dimensionless  sand  transport  distribution  for  the  field 
data  is  compared  with  results  from  the  numerical  model  in 
Fig.  (6.8).  The  measured  maximxm  net  longshore  sand  trans¬ 
port  occurs  on  the  beach  face  at  x  -  20  m.  The  hypothesis  that 
the  predicted  longshore  transport  distribution  is  proportional  to 
the  net  sediment  transport  distribution,  calculated  from 
the  profile  changes  does  not  appear  to  be  the  case.  Signifi¬ 
cant  profile  changes  occur  at  the  beach  face  which  are  not 
included  in  the  model  calculations.  It  is  also  possible 
that  the  sand  eroded  from  the  beach  face  was  moved  in  the 
offshore  direction  outside  100  m  and  was  not  included  in  the 
measured  profiles,  which  would  also  invalidate  the  original 
hypothesis . 

D.  THE  EMPIRICAL  COEFFICIENT  B  AND  LONGSHORE  SAND  TRANSPORT 

s 

The  empirical  coefficient  B  ,  as  described  by  Thornton 

S 

(1973) ,  has  a  magnitude  of  order  (1)  and  is  equal  to 
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Figure  (6.8).  Normalized  sand  transport  distribution 
inferred  from  profile  changes  vs. 
numerical  model  results  with  variable 
depth . 
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(6.1) 


where  e  is  the  transport  efficiency  of  the  available  wave 
power  and  fj^  is  the  dynamic  friction  factor  associated  with 
sediment  transport.  The  factor  f^  is  equal  to: 

^ ,  _  sediment  shear  stress _ 

^h  normal  stress  at  the  bottom  ' 

while  b  is  a  proportionality  coefficient  between  the  ratio  of 
sand  velocities  (^) ^  and  the  square  root  of  the  water  particle 
velocity  ratio 

Because  the  data  do  not  provide  information  about  each 

individual  parameter,  all  parameters  are  combined 

in  B  ,  which  is  determined  empirically  by  model  calibration, 
s 

The  total  longshore  sediment  transport  rate  (Q^)  in 
(m^/day)  is  calculated  by  integrating  the  area  under  the  sand 
transport  distribution 

x*180  _ 

The  calculation  based  on  profile  changes  does  not  give  the 
total  longshore  sand  transport  but  gives  the  net  transport 
which  is  much  less  in  this  case  than  the  longshore  transport 
predicted  by  the  model.  Therefore,  equation  (2.5)  suggested 
by  CERC  (1977)  is  used  as  a  reference  to  calibrate  the  long¬ 
shore  sand  transport  model. 
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The  calculations  of  the  longshore  sand  transport  rate 
require  the  wave  condition  at  the  breaker  line  ,  hj^ 

and  Field  data  give  equal  to  0.9  8  m.  Fig.  (6.2a) 

and  hj^  to  be  2.96  m  knowing  the  beach  slope  besides  the  angle 
cTj^  equals  six  degrees  from  the  refraction  diagram,  Fig.  (4.4). 

calculated  using  equation  (2.6)  gives  4600  m^/day,  which 
requires  to  be  0.18. 

A  summary  of  the  results  if  given  in  Table  III.  The 
difference  between  the  two  cases  of  depth  (constant  and 
varying)  for  the  analytical  model,  was  found  to  be  insignifi¬ 
cant.  A  similar  result  was  found  for  the  numerical  model. 


The  calculated  using  the  numerical  model  is  less  than 

calculated  from  the  analytical  model,  because  the  numerical 

model  used  the  actual  bottom  profile.  The  numerical  and 

the  analytical  models  gave  essentially  the  same  maximum 

cross-shore  transport  distribution  q  ,  for  the  same  case 

^max 

of  either  constant  or  varying  depth.  Fig.  (6.6)  and  Fig. 

(6.7)  . 

Johnson  (1953)  estimated  the  longshore  sand  transport 
at  Santa  Barbara  Beach  to  be  220,000  m^/year,  or  about  600 

3 

m  /day  on  the  average.  The  predicted  by  the  model  gives 
almost  eight  times  the  above  estimated  average  daily  rate. 

The  high  value  of  predicted  by  the  model  reflects  the 
storm  conditions,  which  caused  significant  longshore  transport 
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Summary  of  the  Results 


VII.  SUMMARY  AND  CONCLUSIONS 

The  principle  objective  of  this  study  is  to  predict 
the  longshore  sediment  transport  distribution.  Special 
emphasis  was  placed  on  testing  an  analytical  and  nxamerical 
model  against  field  data  acquired  from  an  experiment  at 
Leadbetter  Beach,  Santa  Barbara  during  the  intensive  storm 
period  17  to  18  February,  1980. 

Field  data  show  a  strong  longshore  current  coupled  with 
significant  changes  in  the  beach  profile  with  no  evidence 
of  on-offshore  sand  movement  (bars  or  troughs)  along  the 
profiles.  Therefore,  it  is  hypothesized  that  the  changes 
in  the  profiles  are  due  to  the  variations  of  sand  transported 
alongshore  by  the  observed  longshore  current,  and  that  the 
on-offshore  sand  transport  is  negligible.  The  net  sand 
transport  calculated  from  the  profile  changes  is  expected 
to  have  a  distribution  similar  to  the  cross-shore  sand 
transport  distribution. 

Initial  values  for  the  models  of  H  _  ,  a  and  f  are 

rms^  o  p 

obtained  from  the  field  data.  Model  coefficients  y  =  0.44, 

B  =  1.1  and  c^  =  0.005  are  determined  by  model  fitting.  The 
developed  model  is  capable  of  predicting  the  increase  in  wave 
height  due  to  shoaling  and  the  decrease  after  wave  breaking. 
The  model  predicts  the  longshore  current  (V)  and  the  longshore 
sand  transport  (q  )  at  each  location.  A  maximum  q  is 

9  S 
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A 


predicted  at  distance  x  =  0.8  which  agrees  with  Komar 
(1977a).  The  experimental  transport  coefficient,  is  found 

to  be  0.18  instead  of  0.08  by  Thornton  (1973).  The  maximum 
longshore  current,  ,  occiirs  at  approximately  the  same 

location  as  (q  )  ,  indicating  that  the  sand  transport 

distribution  is  primarily  determined  by  the  longshore  ciorrent 
distribution. 

The  net  change  in  the  beach  profiles  used  to  infer  the 
net  sand  transport  distribution,  shows  maximum  at  the 
beach  face.  Since  the  sand  transport  model  does  not  consider 
the  swash  region,  the  net  profile  change  distribution  does 
not  agree  well  with  the  predicted  cross-shore  distribution, 
and  the  original  hypothesis  does  not  appear  to  be  valid.  As 
a  result,  sand  transport  mechanics  in  the  swash  zone  appear 
to  be  important  in  determining  why  the  maximum  sand  trans¬ 
port  occurs  at  the  beach  face.  The  swash  is  driven  primarily 
by  infragravity  waves  (long  waves)  due  to  surf  beat  and  edge 
waves  which  have  not  been  included  in  the  surf  zone  models 
where  only  the  depth-limited  sea-swell  frequency  band  of 
waves  are  considered.  Therefore,  further  studies  need  to 
model  the  swash  zone  in  order  to  describe  the  onshore  water 
wave  boundaries  in  a  more  realistic  manner. 
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